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ABSTRACT 

Context. Tidal dissipation in planets and stars is one of the key physical mechanisms driving the evolution of star-planet and planet- 
moon systems. Several signatures of its action are observed in planetary systems thanks to their orbital architecture and the rotational 
state of their components. 

Aims. Tidal dissipation inside the fluid layers of celestial bodies are intrinsically linked to the dynamics and the physical properties 
of the latter. This complex dependence must be characterized. 

Methods. We compute the tidal kinetic energy dissipated by viscous friction and thermal diffusion in a rotating local fluid Cartesian 
section of a star/planet/moon submitted to a periodic tidal forcing. The properties of tidal gravito-inertial waves excited by the pertur¬ 
bation are derived analytically as explicit functions of the tidal frequency and local fluid parameters (i.e. the rotation, the buoyancy 
frequency characterizing the entropy stratification, viscous and thermal diffusivities) for periodic normal modes. 

Results. The sensitivity of the resulting possibly highly resonant dissipation frequency-spectra to a control parameter of the system 
is either important or negligible depending on the position in the regime diagram relevant for planetary and stellar interiors. For 
corresponding asymptotic behaviors of tidal gravito-inertial waves dissipated by viscous friction and thermal diffusion, scaling laws 
for the frequencies, number, width, height and contrast with the non-resonant background of resonances are derived to quantify these 
variations. 

Conclusions. We characterize the strong impact of the internal physics and dynamics of fluid planetary layers and stars on the 
dissipation of tidal kinetic energy in their bulk. We point out the key control parameters that really play a role and demonstrate how it 
is now necessary to develop ab-initio modeling for tidal dissipation in celestial bodies. 
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1. Introduction and context 


Tides have a strong impact on the evolution of star-planet and 
planet-moon systems over secular time scales. Indeed, because 
of the dissipation of the kinetic energy of flows and displace¬ 
ments they induce in celestial bodies, they drive their rotational 
and orbital evolutions (see e.g. |Ogilvie| ( |2014| l and references 
therein for the cases of stars and giant planets and |Correia &] 
Laskar (2003[); Correia et al. ( 2008|l for telluric planets). In the 


case of the Earth-Moon system, their causes and effects are 
now strongly constrained and disentangled thanks to satellite 


altimeter high-precision observations of ocean tides (Egbert 


|& RaypOOOl |2001[ |Ray et aL]|2001[ ). In the Solar System, the 
actions of tides are detected and estimated from high-precision 


geodesic or/and astrometric observations (see e.g. Konopliv & 
Yoder ( I996|l for Venus, Williams et al. (|2014[) for the Moon, 


Lainey et al.|( 2007 1; [Jacobson ( 2010| l; Konopliv et al. ( |2011) l 
for Mars, Lainey et al. (|2009[ 26\2\ for Jupiter and Saturn and 


the attempt by [Emelyanov & Nikonchuk| P013|l for Uranus). 


Einally, a new extremely important astronomical laboratory to 
explore and constrain the physics of tides is constituted by the 
numerous exoplanetary systems discovered over the last twenty 
years ( Mayor & Queloz|1995[ Perryman|201 1 1 . Indeed, they are 
composed of a large diversity of planets (from hot Jupiters to 
super-Earths) and host stars while their orbital architecture and 
the configuration of planetary and stellar spins strongly differ 
from the one observed in our Solar system (e.g. Albrecht et al.| 
|2012[ |Eabrycky et al.|[2012[ Valsecchi & Rasio 2014) 1. In this 


context, the understanding of the tidal formation and evolution 
of planetary systems is one of the most important problems of 
modem dynamical astronomy (e.g. Laskar et al.|2012[|Bolmont| 
et al.|2012] Eerraz-Mello|20 1 3' i while the needed understanding 


and quantitative prediction of tidal dissipation in celestial bodies 
is still a challenge (e.g. Mathis & Remus||2013[ Ogilvie||2014| 
for complete reviews). 


We owe the first theoretical work about a tidally deformed body 
to Lord Kelvin ( |Kelvin|[l863| l. Then, a physical formalism has 
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been elaborated by Love, who introduced the so-called Love 
numbers ( |Love| |1911[ ). The second-order Love number (k 2 ) 
measures the inverse of the ratio between the tidal potential and 
the resulting linear perturbation of the gravitational potential 
at the surface of a body hosting a companion. It quantifies its 
quadripolar hydrostatic elongation along the line of centers. The 
introduction and estimation of tidal dissipation then became 
more and more important with the introduction of the tidal 


quality factor Q for the Solar system (e.g. MacDonald 1964 


|Kaula|[T964[ |Goldreich & Soter||1966| l. This general parameter 
is a simple and useful quantity that allows us to simplify the 
numerical simulations of planetary systems over secular time 
scales. In particular, the constant Q model, the so-called Kaula’s 
model ( |Kaula||T964| l, allows us to take into account internal 
dissipation with one parameter only. The factor Q is applied 
to rocky/icy and fluid bodies and correspond to a simplified 
Maxwell’s model, which assimilates a celestial body hosting 
a companion to a forced damped oscillator. It evaluates the 
ratio between the maximum energy stored in the tidal distortion 
during an orbital period and the energy dissipated by the friction. 
Therefore, the global response to a tidal perturbation depends 
on a mean elasticity, which models the restoring force acting 
on tidal displacements, and a friction, l/Q representing the 
damping coefficient of the system (i.e. the imaginary part of the 
second-order Love number; for a detailed study, see |Greenberg| 
2009 |l. From this approach, other models have been developed 


in order to take into account the possible dependence of Q 
on the main tidal frequency such as the constant time lag 


model (|Darwin| 1 879| [ATexanderl 1 973| |Singer| 1968t |Hut| 1981 [ 
[Leconte et al.|2010|l. 


As the tidal quality factor is inversely proportional to the 
internal dissipation in a body, it is determined by the physical 
properties of this latter. Thus, the internal structure of stars, 
planets and satellites is determinant for the study of tidal effects. 
In this context, recent works have derived the tidal quality factor 
as a function of x and rheological parameters. Examples for 
rocky/icy planetary layers can be found in jEfroimsky & Laineyj 
(|2007|l and Henning et al. ( 2009|l (see also Tobie et al. 2005 [ 
Efroimsky 2012[ Remus et al.||2012bl jCoireia et al.||2014| l. 


Efroimsky & Lainey ([2007 1 illustrate the important role played 


by the internal properties of the perturbed host body in the 
orbital evolution of its system. They focus on the equilibrium 
tide in the case of a circular orbit, which designates the com¬ 
ponent of the strain rotating with the perturber at the frequency 
X = 2(n - Q.) in the frame attached to the central body, n being 
the mean motion of the perturber in a non-rotating frame and Q 
the spin frequency of the primary. The other components corre¬ 
sponding to eigenmodes of oscillation excited by tides constitute 
the dynamical tide. It is particularly important for fluid planetary 
regions and stars in which it can dominate the equilibrium tide. 
It leads to a resonant tidal dissipation that varies over several 
orders of magnitude as a function of their physical properties, 
evolutionary states and x (see Ogilvie 2014 and references 


therein). Applied to an orbital system, a tidal quality factor 
following this resonant behavior makes the orbital parameters 


evolve erratically along secular timescales (Auclair-Desrotour 


jet al.|20f4| ). An interesting example is given by the case of the 
Earth-Moon system. Indeed, most of tidal dissipation on the 
Earth is due to tidal waves in oceans ( [Egbert & Ray|20()0 1 where 
the barotropic tide is converted non-resonantly into internal 
waves because of non-trivial bottom topography (Pedlosky 


1982jl. The corresponding frequency-spectrum of the dissipation 


is then an highly complex function of x while its dependence 


Article number, page 2 of|24| 


on the rotation of the planet allows us to match tidal evolution 
results with geological data. This important result would not be 


possible with constant Q or constant time lag models (Webb 


|1980|[T982) i. During the past five decades, numerous theoretical 
studies have been carried out to characterize and quantify the 


tidal energy dissipated in fluid celestial bodies (see e.g. Zahn 


1966a|b|c| 1975 1977 19891 Ogilvie & Lin 

o 

o 

2005 

Ogilvie & Lin 

2007[ Remus et al. 2012a Ce 

rron et al. 

2012 


stars and the 


envelopes of giant gaseous planets (see |Ogilvie|[2014| ). They 
study global models for tidal waves and dissipative mechanisms 
in fluid regions excited by a tidal perturber that result from com¬ 
plex actions and couplings of rotation, stratification, viscosity 
and thermal diffusion. The kinetic energy of tidal waves is thus 
dissipated by viscous friction and thermal diffusion. Magnetism 
also intervenes in stars and some planets that introduces Ohmic 
diffusion in addition to the two previous dissipative mechanisms. 


Therefore, tidal perturbations take the form of waves in fluid 
regions: acoustic waves driven by compressibility, inertial 
waves driven by the Coriolis acceleration, gravity waves driven 
by the Archimedian force, Alfven waves driven by magnetic 
forces. Given their high frequencies, acoustic waves are only 
weakly excited by low-frequency tidal forcing and shall be 
ignored. Tides will rather excite gravito-inertial waves. These 
mixed waves combine the second and third families enumerated 
above. Inertial waves are caused by rotation. The Coriolis 
acceleration acts as a restoring force and their frequencies are 
smaller than the inertial one 2Q. Gravity waves propagate in 
stably-stratified fluid regions, with frequencies bounded by N 
(the Brunt-Vaisala frequency) that depends on the gradients of 
the specific entropy. Their restoring force is buoyancy. Alfven 
waves result from the presence of a magnetic field and propa¬ 
gate in stars and planets, gravito-inertial waves thus becoming 
magneto-gravito-inertial waves. Erom now on, magnetism is 
neglected in a first step and they are thus not studied in this work. 


Given the complex dynamical resonant tidal response of fluid 
planetary layers and stars described above and its important as- 
trophysical consequences, it should be understood and charac¬ 
terized systematically. To reach this goal, two approches can 
be adopted. Eirst, global realistic models constitute an efficient 
tool to explore the physics of dissipation and to quantify it 
(e.g. Ogilvie & Lin 2004 2007| l. Simultaneously, local sim¬ 
plified models allow to understand in detail complex physical 
mechanisms in action and to explore the large domain of pos¬ 
sible parameters in astrophysics (e.g. Ogilvie 2005 1 Jouve & 
|Ogilvie|2014[ [Barker & Lithwick||2013| [2014 1 . In this frame- 
work, [Ogilvie & Lin[ ( 2004| l proposed in an appendix of their 
work devoted to quantify global tidal dissipation in spherical 
shells for planetary interiors a reduced model based on a Carte¬ 
sian fluid box. The objective of this simplified model was to un¬ 
derstand the behavior of tidal dissipation caused by the (turbu¬ 
lent) viscous friction acting on inertial waves in planetary and 
stellar convective regions. Similarly, [Gerkema & Shrir^ ( [2005a| l 
explored the complex physics of gravito-inertial waves (see also 
[Mathis et al. [20145) 1. Eollowing these previous studies, this work 
develops a new local Cartesian model that generalizes the one 
by Ogilvie & Lin ( [2004| l by taking into account the relative in¬ 
clination between gravity and rotation, a possible stable entropy 
stratification or a super-adiabaticity established by convection 
and viscous and thermal diffusions. By studying this set-up, we 
answer to the following questions: 
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• how does the dissipation due to the viscous friction and heat 
thermal diffusion depends on the fluid parameters (i.e. rota¬ 
tion, stratification, viscosity and thermal diffusivity)? 

• how do the properties of the corresponding dissipation fre¬ 
quency spectrum vary as a function of the latter? 

• what are the asymptotical behaviors relevant for planetary 
and stellar interiors? 

In section 2, the local Cartesian model is presented. Dynami¬ 
cal equations are derived and solved analytically. It allows us to 
get the energy dissipated per unit mass in the box over a rotation 
period. Physical quantities are written as a series of resonant nor¬ 
mal modes corresponding to the harmonics of the tidal forcing 
allowing us to study the spectral response of the box. In section 
3, it is studied on the one hand to identify the four asymptotic 
regimes of tidal waves. On the other hand, we characterize the 
corresponding resonant dissipation frequency spectrum, i.e. the 
positions, width, heights, number of peaks, the level of the non¬ 
resonant background and the sharpness ratio (defined as the ratio 
between the height of the main peak and the level of the non¬ 
resonant background). Each of these properties is expressed as 
a function of the fluid parameters and forcing frequency. These 
scaling laws give a global overview of the properties of the dis¬ 
sipation. In section 4, we examine the particular case of super- 
adiabaticity in convective regions for which the square buoyancy 
(Brunt-Vaisala) frequency is negative. In section 5, we discuss 
our results. Finally, in section 6, we give our conclusions and 
highlight new questions for further studies. 


2. Forced waves in stars and fluid planetary layers 

2.1. The local model to be studied 


Ogilvie & Lin (2004 1 proposed a local Cartesian model describ¬ 
ing the behavior of a fluid box submitted to tidal perturbations 
for convective regions where inertial waves propagate. This 
approach is of great interest because it yields an analytical 
expression for the viscous dissipation of the tidal kinetic energy, 
which is expressed as a function of the fluid properties and 
particularly of the Ekman number, A^Ek = v/ Q being 

the angular velocity of the body, v the local kinematic viscosity 
of the fluid, and L a characteristic length of the convective 
layer. So, it allows us to understand how dissipation depends on 
these physical parameters. Here, in addition to inertial waves 
driven by the Coriolis acceleration, we also take into account 
the possibility of a stable stratification that introduces buoyancy 
as an additional restoring force. This allows us to generalize the 
previous model to study the more general case of gravito-inertial 
waves (e.g. Gerkema & Shrir^|2005a| l. Moreover, we also take 
into account thermal diffusion of heat in addition to viscous 
friction to dissipate tidal kinetic energy. 


We consider a local fluid Cartesian section belonging to a planet 
or a star tidally perturbed at a frequency x'- ^ box of length L, 
such that L <s: R, where R is the radius of the body. The fluid is 
newtonian, of density p, kinematic viscosity v and thermal dif¬ 
fusivity K. The box vertical direction (and gravity g) is inclined 
relatively to the spin vector of the body with an angle 0. The 
angular velocity of the fluid (Q) and the local gravity g are as¬ 
sumed to be uniform and constant. The center of the body and 
of the box are denoted O and M respectively. We use two ref¬ 
erence frames (Fig. [TJ. The global one • {<5, Xe, Ye, Ze) 

rotates with the body and its natural spherical associated unit- 
vector basis is denoted eg, e^^ In this frame, Q = QZe and 


the coordinates of M are (r, 0, (p) in spherical geometry. Then, we 
introduce the frame fixed to the fluid section, R : |m, e^, e^, e^j 
whose unit-vectors are = e^, = -eg, 6;^ = and g = -ge^.. 

Stratification implies a new frequency, specific to gravito-inertial 
waves, the Brunt-Vaisala frequency denoted N and defined by 


N^^-g 


dlogp 

dz 


I d log P 
y dz 


( 1 ) 


where j - (dlnPjdlnp)^ is the adiabatic exponent {S being 
the specific macroscopic entropy), and P and p the radial distri¬ 
butions of pressure and density of the background respectively. 
These equilibrium quantities are supposed to vary smoothly with 
the radial coordinate z compared to the perturbation. For this rea¬ 
son, the Brunt-Vaisala frequency is taken as a constant parame¬ 
ter. Assuming the hydrostatic equilibrium, it becomes 




g\dp ^ gp^ 
p[dz yP 


( 2 ) 


The regions studied are stably stratified (A^ > 0) or convective 
(A^ 0 or A^ < 0). 


Q 



Fig. 1. The fluid box, its reference frame and its position in the planet 
relatively to the spin axis. 


Inertial waves 

Convective Layers 

^ (assuming N^-O) 

- i -1- 1 - 1 - 

tnf.ssL 2Q A/gsL ^«up;ssL 3 tai 3 |y Stratified Layers 

Gravito-inertial waves 

Fig. 2. Gravito-inertial waves frequency range in Convective and in 
Stably Stratified Layers (Ac = 0 and Asll > 0 are the corresponding 
buoyancy frequencies). The expressions of the boundaries ;^'mf;ssL and 
Asiip;SSL are explicitly given in Eq. ( (^ . 


2.2. Dynamical equations 

Submitted to a tidal forcing massic force F = (^Fx, Fy, F^, where 
X = (x, y,z) are the local coordinates of space, the fluid moves in 
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the box periodically and dissipates energy through viscous fric¬ 
tion and thermal diffusion. To obtain the equations that govern 
the dynamics of tidal gravito-inertial waves, we write the lin¬ 
earized equations of motiorQof the stratified fluid on the non- 
traditional /-plane assuming the Boussin escp] and the Co wling 
approximations (Gerkema & Shrira 2005 a[ [Cowling| 1 94T i. First, 
we introduce the wave velocity field u = (m, v, w), where u, v and 
w are its components in the local azimuthal, latitudinal, and ver¬ 
tical directions. Next, we define the fluid buoyancy 


the energy dissipated by thermal diffusion 


n _ D _ P (x, f) 

B — ~ 9 ’ 

P 


(3) 


where p and p are the density fluctuation and the reference back¬ 
ground density, and t is the time. The three linearised compo¬ 
nents of the momentum equation are given by 

dfU) -2Q.s,m6u - — d,p +B + vV^w + F^ 

d,v + 2D.cos6u - — dup +vV^v + Fu , (4) 

^ 1,2 

dtU - 2Qcos 6v + 2Q. sin6u; = — d^p -i- vV m -i- F^ 

P 

where p is the pressure fluctuation and + dy^y + 

Next, we write the continuity equation in the Boussinesq approx¬ 
imation 

V ■ u = d^w + dyV + dxU = 0. (5) 

Finally, we get the equation for energy conservation 

d,B + N^w = kV^B. (6) 

By taking the scalar product of the momentum Muation (Eq. 
with u and multiplying the heat equation (Eq. by BjN^, we 
obtain the energy equation for the box 


d, (^Ec + Epj^- r V ■ (p'u) dV 

«y V 

where we have introduced 
• the kinetic energy 


H- -H D' 


.them 


^forcing^ 


1 

2 


I pu^dV, 

Jv 


the potential energy associated with stratification 


\e,. = \ 


X 


if 

.f2 _ I 


[Ep = 0 if = 0 


the energy dissipated by viscous friction 
J p(vu-V2u)dy, 


( 8 ) 


(9) 


( 10 ) 
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= j p^j^BV^BjdV if 

Utherm^Q jf ^2 ^ q 


(11) 


( 12 ) 


• the energy injected by the forcing 

^forcing ^ r p (u . F) dV, 

Jv 

where V is the volume of the box. The first term on the right- 
hand side of Eq. (j7]i is the divergence of the acoustic flux. 

We introduce dimensionless time and space coordinates and tidal 
frequency, normalized buoyancy and force per unit mass 


T - 2Llt, X - r=—, Z - 0 ) - —, 

L L L 2Q 

b=— f= — 

2Q’ 2G' 

The linearized Navier-Stockes equation becomes 

where we recognize the Ekman number 

V 

Nek = 


-Vp - AfikV u - b = f. 


(13) 


(14) 


(15) 


2GL2’ 

which compare the strength of the viscous force to the Coriolis 
acceleration. 


The equation of heat is written 
drb + Aw - NdiffV^fe, 


(16) 


where the right-hand side corresponds to thermal diffusion. We 
introduce here two new control parameters 


A = (A) and = 


* In this work, we do not take into account the nonlinear interactions 
between tidal waves while those with turbulent flows such as convection 
are treated as an effective viscosity ( [Goldreich & Keeley||1977{ |Zahn| 
1989||Ogilvie & Lesur|2012|l. 

^ lidal frequencies are most of the time small compared to the Lamb’s 
frequency. This allows us to filter out acoustic waves and to adopt the 
anelastic approximation in which V • (p u) = 0. In the case where the 
background density p is uniform and constant, it reduces to the Boussi¬ 
nesq approximation in which V • u = 0. 


(17) 

2QL2’ 

A is the square ratio of the characteristic frequencies of the sys¬ 
tem (2Q and N). This parameter is bound to the nature of tidal 
waves (e.g. [Gerkema & Shrir^ 2005 a|b| see also jBerthomieuj 
jet al.|(|l978]l andjProvost et al.|(|1981|l): 

• A < 0 (i.e. < 0) corresponds to inertial waves, 

• 0 < A < 1, corresponds to quasi-inertial waves, 

• A > 1, corresponds to gravito-inertial waves. 

The parameter Adiff compares the strength of the Coriolis and 
thermal diffusion terms by analogy with the Ekman number de¬ 
fined above (Eq. [T5] l. Einally, we introduce the Prandtl number 
of the fluid Pr that characterizes the relative strength of viscous 
friction and thermal diffusion: 

V Apk 

F, = - = (18) 

K Adiff 

Waves are damped by viscous diffusion if Fr > 1; otherwise, 
thermal diffusion dominates. 

Erom now on, we assume that the unknown quantities vary 
with X and z only {x e [0,L], z e [0, L]). This assumption has 
no impact on qualitative results and allows us to compare our 
calculations to the ones made by Ogilvie & Lin ( 2004|l directly. 
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2.3. Velocity field and dissipation 

We now solve our linearized first-order system (Eqs. [I4l[5land 
m to get the velocity field and the energy dissipated by viscous 
friction and thermal diffusion. The tidal perturbation being peri¬ 
odic in time, we expand the quantities in Fourier series: 




, Uy = %[v{X,Z)e^‘‘^^ 


u;(A,Z)e-'"^ 

, p = !R[i/^(A,Z)e-'"^' 


/(A,Z)e--^ 

, fy^%\g(X,Z)e-“^^ 



, b ^ %[b{X,Z)e-‘p 


The spatial functions are themselves expanded in periodic spatial 
Fourier series in X and Z as in|Ogilvie & Lin|(f2004|l: 


u — 'y 

h 'y 


i2K{mX+nZ) 


ilnimX+nZ) 






i2Tr(mX+nZ) 


ie 


,i2n{mX+nZ) 


i2n(mX+nZ) 


i27T{mX+nZ) 


g^'^gn. 


( 20 ) 


i2nimX-^nZ) 


i27T{mX+nZ) 


For the sake of simplicity, the horizontal and vertical wave 
numbers {m,n) e Z*^ are not written under summations. Note 
that by choosing periodic boundary conditions, we consider 
the volumetric excitation of normal modes ( |Wu|2005j |Braviner| 
& Ogilvie 2015|l. The case of periodic wave attractors can be 


treated by introducing reflexions on boundaries of the box (e.g. 
|Jouve & Ogilvie|2014|l. The previous system becomes: 


-itOUmn - COS 6v,„„ + siu OWmn 

-iuVmn + COS Ou^n 

sin OiifYifi 


—imAi//„j„ — E {in^ + 
'^fmn 

-E {m^ + Vmn + g,„„ 
—inAi//mn - E {irP' + 


mUrnn + nWmn - 0 

-iL0b,„„ + AWmn = K {m^ + rP’') bmn 

( 21 ) 

which is parametrized by the frequency ratio A, the colatitude 0 
and the numbers: 


^ 27r^v ^ ItPk 


and 


A = 


n 


( 22 ) 


E and K are both dimensionless diffusivities. E is proportional 
to the Fkman number Nsk (Eq. [TSl l of the fluid and K to its ana¬ 
logue Adiff (Fq. for thermal diffusion. Similarly, we get for 
the Prandtl number (Fq. [TSll 


Pr = 


E 

K' 


(23) 


Finally, A weights the pressure variations. This last parameter 
does not intervene in the expressions of the velocity field and of 


the perturbation of buoyancy (Eqs. 29 and 131). Therefore, both 
viscous dissipation and thermal diffusion do not depend on it. 
The equations yield two complex dimensionless frequencies, be¬ 
cause of viscous and heat diffusion, that are functions of E and 
K: 


(b - oj + iE (irP + rP^ and cb - oj + iK {m^ + rP^ . (24) 

Initially, assuming that f = 0, we get the dispersion relation of 
the viscously and thermally damped gravito-inertial modes: 

rPcoP'9 ifpA (b 

d) - —r- Y -T—. (25) 

For slightly damped modes (£ <; 1 and K ^ 1), we identify 
in the second member the wave number k = (kH,Q,kv), with 
kff - mjE and ky - njL (see |Gerkema & Shrira|2005a) . Indeed, 
neglecting E and K, it simplifies to 


irPN^ + tP (2Q. cos 6) 


that can be written 


(26) 




(27) 


where we recognize the dispersion relation of gravito-inertial 
waves with the respective contributions of inertial and gravity 
waves. 

By letting m and n tend to infinity successively, we identify the 
approximative boundaries of the frequency range, Xin! and Xi 
When 2Q < N (i.e. A > 1), A'inf ~ 2fi and;^fsup ~ N (see Fig. jS] 
The exact expressions of xm and;^fsup are given by|Gerkema 
[Shnra| ( |2005a) : 


Tint 


Tsup 



/12 - (2o-Nf 

1/2 

\ 2 

A 

A2 - (2crNf 

1/2 


2 


(28) 


where A - + cr^ + cr^. Here, cr = 2Gsin0 and ctj = 

2Q cos 6 sin a, a = 0 corresponding to the propagation direc¬ 
tion of gravito-inertial waves in the horizontal plane along the x 
axis (see the bottom-left panel of Fig. [^l. At the end, we obtain 
the coefficients of the velocity field u, 


^mn 


i(b {nfmn - mhfnn) " ^ COS Og^n 

{m? + rP-) ap- — rP cos^ 6 — AtrP — 

OJ 


n cos 6 (nf„„ - mh^„) + i 

/ 2 2\ ~ 2lm2 

[m + n ]u - 

0mn 


^ ^ OJ 



{nP + rP) aP — fP cos^ 6 — AnP — 

cb 


i(x) (jlfmn ^ COS 

—m -, 

{mP + rP) dP — rP cos^ 6 — AnP — 

u 
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(29) and its viscous and thermal contributions. 


of the pressure p. 




1 

A 


(w/„„ + i cos 6g 

mn) 


n sin 0 + im — 




(m^ + n^) (xP- - rP- cos^ 9 - Am^ — 


hmn sin OCj + in {cP - cos^ 


U) 


(trP + P) dP — cos^ 9 — ArP — 
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Fourier series expansions given in Eq. (20 1 , this expression be¬ 
comes: 
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Using once again Eq. ([20|), it becomes: 
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The total dissipated power is the sum of both contributions: 
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that will be studied instead of the power (D). The energy (^ of 
Eq. 


Eq. 


221 . 


37 1 only depends on the control parameters A, E and K (see 
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Note that all the coefficients in Eqs. ( [29| , (301 and ( [3T] l have the 
same denominator. It corresponds to the left-hand side of the 
forced momentum equation (Eq. [I4j l. The tidal excitation only 
affects their numerators. 

As pointed out before, the viscous dissipation of gravito-inertial 
waves over a tidal period is deduced from the velocity 

field (Eq.[T0]l. It is provided by the quadrature of the local mean 
dissipation in the whole box. Using our dimensionless coordi¬ 
nates, we define the viscous dissipation per unit mass: 


2.4. Spectral response 

We have computed analytical formulae of the energies dissipated 
by viscous friction and thermal diffusion. Considering the dis¬ 
persion relation Eq. ( [25| ) and the resonances involved, we know 
that ^ significantly varies with the tidal frequency o). As shown 
by Eq. ( [38] l, its behavior is determined by the position 9 of the 
box and the dimensionless numbers A, E and K. To be able to 
make our results comparable with the model in the appendix of 
[Ogilvie & Lin|(|2004|), we choose the same academic forcing: 


fmn ~ All)’’ — 0 and bffjfi — 0. 

4 \m\ 


(40) 
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where = (1/T"^)V^ and (■■■) is the average in time. Using 


Y +n^){Kn\+Kn\+ymn\)’ ( 33 ) 


as demonstrated by |Ogilvie & Lin| ( |2004| l. Similarly, the power 
dissipated by heat diffusion over a tidal period can be derived 
from the buoyancy perturbation (Eq.fTT). Using our dimension¬ 
less coordinates, we define the heat diffusion per unit mass: 
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Einally, we introduce the energy dissipated over a rotation period 
of the planet, denoted 


(37) 


We verified that the forms of the {fmn,gmn,hmn] coefficients do 
not affect properties of the spectra that are governed by the 
left-hand side of the forced momentum equation (Eq. [T4| , i.e. 
the position and the width of the resonances in the frequency- 
dissipation spectra. However, the spectral dependence on m and 
n modulate their number, their height, and the amplitude of the 
non-resonant background. Nevertheless, the contrast between 
the latter and the height of the main resonance, which will be 
called the sharpness ratio and evaluates the relative amplitude 
between dynamical and equilibrium tides, are less affected. 
Moreover, we verified that the different main asymptotic behav¬ 
iors of the system remain unchanged. In future works, it would 
be important to implement a more realistic treatment of the 
forcing term, as done for example by [Barker & Lithwick ( |2013| l 
for local numerical simulations of the tidal elliptic instability. 

In Eigs. and we plot the frequency-spectra of for 
various sets of parameters. The abscissa measures the dimen¬ 
sionless frequency w = ^nd the vertical axis the local 

viscous dissipation per unit mass (in logarithmic scale). Figure]^ 
corresponds to pure inertial waves (A = 0) in a box located at 
the pole (9 = 0). The plots illustrate most of the results of this 
section. First, the cutoff frequency for inertial modes, tOc - I 
(see the dispersion relation in Eq. [25]l, clearly appears. Beyond 
w = 1, the dissipation decreases by several orders of magnitude 
and corresponds to the non-resonant background. Second, we 
can observe the sharp dependence of the dissipation to the 
tidal frequency x- ^ varies with the Ekman number E. The 
spectrum is regular for high values of E, that correspond to 
waves damped by a strong viscous diffusion. At the opposite, 
the number of resonances increases when E decreases. 

Eigure gives an overview of the dependence of and 

^therm contTol parameters. Each plot corresponds to a 

particular regime of tidal dissipation. The number of peaks 
is obviously correlated to the Ekman number and we retrieve 
the spectra of Eig in the top-left and top-right plots where 
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dissipation by viscous friction predominates. The ratio of the 
contributions, ^''iscy'^therm^ depends on the Prandtl number, the 
dissipation being mainly due to viscous friction for high values 
of Pr and to heat diffusion for low ones. The cases of the 
bottom-left and bottom-right plots show that the critical Prandtl 
number, where both contributions are comparable, varies with A. 

After this first qualitative approach, it becomes necessary to de¬ 
velop a quantitative physical description of the spectra taking 
into account both dissipation mechanisms. In a previous work, 
we had pointed out the importance of the role played by reso¬ 
nances in the evolution of planetary systems. Indeed, the varia¬ 
tions of the distance between a star or a fluid planet and its planet 
or satellite are directly linked to the height and width of the peaks 
(see |Auclair-Desrotour et al.|2014| l. It is now possible to under¬ 
stand how the latter depend on the physical fluid parameters. 


3. Properties of the resonances 

Numerical results show that the fluid behaves like a bandpass fil¬ 
ter with clearly identified cutoff frequencies, cuinf and Wsup, com¬ 
mented further (see Eq. [28l l. A typical dissipation spectrum is a 
batch of resonances presenting a large range of positions, widths 
and heights, and an apparent non-resonant background. Thanks 
to our analytical formulae of the viscous and thermal dissipations 
(Eqs.[^and|^l, these properties can be established explicitly 
as functions of the colatitude 6, the ratio of squared frequencies 
A, the Ekman number E and its analogue for thermal diffusion K 
(and the corresponding Prandtl number Pr). Thus, we will treat 
hrst the eigenfrequencies Umn and the widths that do not de¬ 
pend on the forcing, and in a second step the heights and 
the non-resonant background i/bg- At the end, we will look at 
the number of resonances Akc and the sharpness ratio E, which 
is the ratio between the height of the main peak and the mean 
non-resonant background inside the frequency range [cuinf, Wsup]- 


3. 1. Position and popuiation of resonances 

The spectra of Eig. reveal a well-organized structure. Each 
peak is surrounded by two others of the nearest higher orders. 
The fractal pattern, accentuated by the particular forcing chosen, 
is apparent in the bottom-right plot (E - 10^^). In this subsec¬ 
tion, we retrieve this singular structure analytically, by studying 
the expression of ^ (Eq.[38|). The energy dissipated is written as 
an infinite sum of terms associated to the degrees m and n. 
As (u>) (tu) (tu) i^o), 171 and ii being 

natural integers, four terms of this series dominate all the others 
near a resonance. So, we assume: 


minima of the denominator dm„. Then, we factorize the expres¬ 
sion given in Eq. ( |43] l by 1/d), which varies smoothly like the 
numerator. It yields the polynomial P: 

P (w) - 0 )^ + - 2a) + (a^ - 2/ly) 0 / +y^, (44) 

where 

a 

' J3 

r 

P (ai) have the same minima as d,„n. So, by differentiating P (w), 
we get: 


a, p and y are positive real coefficients: 

/ . -N 2 rP cos^ 6 + irPA 
= E(E + 2K)(m^ + nA + -^^- 

= (K + 2E) {nP + rP) 

= KE^ {nP + rP'f + rP cop- OK + npAE 


(45) 


P'(a))^a)Q(oj), (46) 

Q being an even fourth degree polynomial. Thus, solutions sat¬ 
isfy the equation: 


e(tu) = 0, (47) 

and exist only if the coefficients a, p and y satisfy the condition: 

+6Py+pP - A-aP^ > 0. (48) 

Given the form of the polynomials P and Q, which involve only 
even terms, eigenfrequencies are symmetrical with respect to the 
vertical axis (ai = 0). Therefore, the domain studied can be re¬ 
stricted to the positive interval [0; -i-oo[. There is a single physical 
solution if 2a > p^ or a^ < 2py, 


^mn ~ 


1 


2a - P'^ + ■Jp‘^ + a^ - 4aP^ + 6Py 


(49) 


which corresponds to a minimum of P, so a maximum of We 
have pointed out that and have the same denominator. 
Therefore, they also have the same eigenfrequencies. In a first 
time, consider the bi-parameter case of Ogilvie & Lin (2004i: 
A - 0, K - 0 (the solution only depends on the position 9 and 
the Ekman number E). The eigenfrequency tOmn becomes: 


^v.sc ^ 4^v.c ^ ^ „2) ^ |^2^J ^ 1^2 ^ ( 41 ) 

^therm ^ 4^*erm ^ 8^^-2£p-l ^^2 ^ „2^ _ (42) 

and have the same denominator. 
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Considering that the numerators varies smoothly compared to 
d,„„, we deduce the eigenfrequencies a),„„ (with (m,n) e N*^) 
from the dispersion relation (Eq. [25|l. They correspond to the 


cos^ 9 


■ E"^ {nP + rP)^. 


(50) 


Beyond a critical rank, there is no minimum any more. The first 
harmonics only are resonant because they verify the condition: 


(jfp + rP) 



(51) 


Thus, there is obviously a maximum number of resonances that 
decreases with E and 9 in the bi-parameter case. However, we 
will see that this maximum does not correspond to the effective 
number of resonances. This latter is actually strongly constrained 
by the level of the non-resonant background which is the cause 
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Fig. 3. Frequency spectra of the energy per mass unit locally dissipated by viscous friction for inertial waves (A = 0) at the position 9 = 0 
for K = 0 and different Ekman numbers. In abscissa, the nor malized frequency a) = ^/20. Top left: E = 10“^. Top right: E = 10 Bottom left: 
E = 10“'*. Bottom right: E = 10“^. The results obtained by Ogilvie & Lin j2Q04| are recovered. Note that the non vanishing viscous dissipation 
at tu = 0 is the one of the geostrophic equilibrium. 


of the variations observed in Fig s.[^ and|4] By plotting the analyt¬ 
ical eigenfrequencies Umn (Fig. 1^ we retrieve the global struc¬ 
ture of spectra. Each blue point corresponds to a mode (m, n). 
The characteristic rank of the harmonic k = maxjlml, \n\] is in¬ 
troduced to reduce the doublet (m, n) to one index only. It repre¬ 
sents an approximative degree of the resonance. The cutoff fre¬ 
quency of inertial waves Uc - cos 6 distinctly appears on Fig. 
1^ given that E ^ 1. It can also be deduced from the simplified 
expression of oJmn'- 


n 

Umn COS 6. (52) 

yirE -H rP- 


Thus, for the main resonances (k < 10), the coefficients a, ji and 
y become: 

rP coP 6 + m^A 
a = -r--- 

' p = (K + 2E){in^ + rP) 

y = rP cos^ 6K + rrpAE 
and the frequency of the resonance (m, n), 

= Va(l -He), (54) 

where e is a second order bias: 


The plots representing the eigenfrequencies (Fig. expose 
the symmetries of the dispersion relation in the quasi-adiabatic 
approximation, that we develop from now on. 


The quasi-adiabatic approximation means that the diffusivities 
E and K are first-order infinitesimals with respect to cos 6 and 
Va. Such an assumption corresponds to low viscous and thermal 
diffusivities (v and k) relevant in stellar and planetary interiors. 



(55) 


Therefore, when E and K are small and can be neglected, reso¬ 
nances are located at the eigenfrequencies 


cos^ 0 + AirP 


(56) 
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Fig. 4. Frequency spectrum of the energy per mass unit dissipated in different regimes (f) (continuous red line), and of its viscous and thermal 
contributions (dashed green line) and ((joshed blue line). The latitude is fixed at 0 = 0. In abscissa, the normalized frequency oj = ;t'/2n. 
Top left: A = 1.0x10 -;£= 1.0x10 ^ A" = 1.0 x 10 ^ Top right: A = 1.0xl0^£= 1.0x10 1.0 x 10 ■*. Bottom left: A = 1.0x10 ^; 

£ = 1.0 X 10 ^ £ = 1.0 X lO T Bottom right: A = 1.0 x 10^; £ = 1.0 x 10 £ = 1.0 x lO T 


that correspond to the dispersion relation of gravito-inertial 
waves in the adiabatic case given in Eq. ( p^ . They give us the 
boundaries of the frequency range in the local model. 


Winf ^ — and Wsup ~ cos ff if N <2Q. (57) 

or the contrary if > 2Q. Note that the exact expressions of 
coiaf and tUsup (Eq. [28| ) are actually a bit more complex than the 
previous ones. In Eigs.j^andj^ the interval [wmf, Wsup] delimits 
peaks and blue points zones. As the impact of E and K is of 
second order, the structure of the batch of resonances only 
depends on A and the colatitude 6 under the quasi-adiabatic 
approximation. This is the reason why the plots of Eig.[^ where 
[E,K) <si jVA, cosflj have the same appearance. By taking 
A = 0, we recover the expression that corresponds to pure 
inertial waves (see Eq.[50|l. 

Moreover, Eq. ( |56] l points out an hyper-resonant case related to 
quasi-inertial waves (A < 1) and characterized by the equality; 


A - cos^ 0. 


When it is verified, all the peaks are superposed at the eigenfre- 
quency: 


ojp - iOc - COS 6, (59) 


and form a huge single peak. So, if A < 2Q (quasi-inertial 
waves), there is a critical colatitude 6o = arccos (N/ (2Q)) at 
which the dissipation spectrum of tidal gravito-inertial waves 
reduces to a single resonance. The frequency range decreases 
when the latitude comes closer to 0q. As there is no such critical 
colatitude for A > 1, usual gravito-inertial waves are not affected 
by this effect. Their frequency range [winf, cUsup] is only larger at 
the equator than at the pole. 

The number of resonances constituting a spectrum can be es¬ 
timated from the expressions of the eigenfrequencies w,„„ (Eq. 
[56] l. Indeed, due to the dehnition of k, there are 2k - I points 
per line on Eig. But the symmetries of Wmn make that some of 
these points have the same positions. Indeed, note that atm'n' = 
if m' - pm and n' - pn with p eZ*. This explains why res¬ 
onances are not as numerous as modes. We denote kc the max¬ 
imal rank of the harmonics dominating the background and 
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Fig. 5. Structures of the frequency spectra of dissipation for inertial waves and gravito-inertial waves, with various sets of parameters. The 
positions of resonances (in abscissa, the normalized frequency w = ^/2f2) are indicated by blue points as functions of the characteristic rank 
k of the harmonics (ordinates). Top left: {B,A,E,K) = (o, 0,10“^, oj. Top right: (0,A,E,K) = (;r/6, 0,10“'',o). Bottom left: (9,A,E,K) = 
[n/6,25, 10-^0). Bottom right: (e,A,E, K) = {n/6,25,0, 10^^). 


the number of effective resonant peaks. The layer of harmonics 
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Nkc 

1 
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7 

11 

19 

23 

35 

43 

55 

65 


Table 1. Numerical comparison between the number of peaks and 
Pk — ~ ^ ~ P' ■ the number of modes k^, k^ being the rank of the highest harmonics, for 

; I (i/iriENMo II jjjg resonances (1 < k < 10). 

k; onme number) ^ ' 


So, the effective number of resonances can be computed using 
the following recurrence series: 


kc 

Nkc^'^Pk- (61) 

c=i 

The values of Nkc and for the first orders are given in Table 
and plotted in Fig. We notice a slight difference between the 
two curves, the effective number of peaks growing slower than 
the number of modes. As a first approximation, we can consider 
that Nkc k^. This assumption will be used to compute the scal¬ 
ing laws of Akc- 

3.2. Width of resonances 

The width of a resonance is a characteristic of great interest. In¬ 
deed, in a previous work (|Auclair-Desrotour et al.|[2014|l, we 


studied a coplanar two-body system of semi-major axis a and 
showed that a was submitted to rapid changes due to resonances 
(see also |Witte & Savonije|1999]|2001||20()2| l. In this work, the 
system is constituted by a fluid rotating planet and a pointlike 
satellite. The variation Aa of the semi-major axis a is related to 
the width at mid-height Ip and sharpness ratio Ep - HplH^^ of 
the p-th resonance encountered. Indeed, Aa/a oc tpEp^, i.e. the 
amplitude Aa linearly varies as a function of the width of the 
peak that causes it, this width being the signature of the dissi¬ 
pation. In this section, we compute an analytical formula of the 
width at mid-height l„„ as a function of the control parameters 
of the system. Similarly to the eigenfrequencies, l^n are fully 
determined by the left-hand side of Eq. ( [I4jl. We suppose that 
the numerator of in Eqs. (pT) or (|42|l, varies smoothly 
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Fig. 6. The real number of resonances (continuous red line) and 
its first order approximation (dashed green line) as a function of the 
rank of the highest harmonics for the main resonances (1 < A; < 10). 


compared to its denominator d„,„ (see Eq. |43] l. Then, the width 
at mid-height is computed from the relation; 


1 


^mn T 2 y — 2 

that can also be expressed: 

^mn 1 ^mn 


r, I ^mn 

P I ^mn T 


2 P (oJ/jin) 


Therefore, we solve the equation: 


(62) 


(63) 


+= 2P(w^„). (64) 

So, the lf„„ are the same for and ^*erm^ eigenfre- 

quencies. In the context of the quasi-adiabatic approximation, 
we obtain: 


Am^K + {in^ cos^ 6 + E 

rP cos^ 6 + AifP 

Looking at the form of this expression, we introduce two 
numbers proper to the mode (m, n), 


Imn — 


{m^ + 


(65) 

critical 


A,„n iS) = ^ COS^ 0 and Pr,mn (0,A) = --——, (66) 

nP A+ A„„ {6) 

that determine asymptotical behaviors. The expression of l^n 
thus becomes; 


7 (2 2 \ ^ (■'^ Afjin) {Pr d" Pr\mn) 

Imn = [rrP + n^j - 

Pr 

A <si A„in characterizes inertial waves and A » A„,„ gravito- 
inertial waves. In the same way, if Py <s: Py-mn, the resonance is 
damped by thermal diffusion and if Py » Py^mn, it is damped by 
viscous friction. We thus identify four possible regimes, summa¬ 
rized in Fig. |7] 


A -H —At, 


1. A A„,„ and Py » Pr,mn' inertial waves damped by viscous 
diffusion; 

2. A » Amn and Py » Py^mn'- gravity waves damped by viscous 
diffusion; 

3. A « A„,„ and Py « Py^ynn'- inertial waves damped by thermal 
diffusion; 

4. A » A,„„ and Py <si Py^mn- gravity waves damped by thermal 
diffusion. 



- 4 - 3 - 2-101234 

log A 

Fig. 7. Asymptotic domains. The areas on the left (a and c) correspond 
to inertial waves, the ones on the right (b and d) to gravito-inertial 
waves. Resonances are damped by viscosity in areas a and b, and by 
thermal diffusivity in areas c and d. The corresponding domains of pa¬ 
rameters are: a: A « A„„ and Py » Py-mn', b: A » A„„ and Py » Py-mn', 
c: A « A„„ and Py « d: A » A„„ and Py « Py-„,„. 

The scaling laws of Tab. are directly deduced from Eq. ( |65] l. 
They illustrate the differences between the four identified 
regimes. First, we focus on inertial waves. If the viscous term 
overpowers the term of heat diffusion, then we are in the case 
studied by |Ogilvie & Lin| p004| l. The width at mid-height 
of resonances linearly varies with the Ekman number. For 
E - 10“^, peaks are larger than for E - 10“^ (see Fig. |^. 
Elsewhere, the width is proportional to K: for a given A > 0, 
the resonances would widen with K as they do with E in the 
previous case. Now, if we examine gravito-inertial waves, they 
behave similarly as inertial waves, linearly widening with E and 
K in the corresponding regimes defined above. Finally, note that 
Imn always depends on a single parameter except in the case 
of inertial waves damped by thermal diffusion, for which the 
square frequencies ratio (A) has also a linear impact. 

Assuming that the resonances all have the same qualitative 
behavior, we plot the width at mid-height /n of the main one 
from Eq. ( |65| ). We thus retrieve in the plots the regimes predicted 
analytically as illustrated in Fig.[^ in particular transition zones. 
Concerning the nature of the damping, the transition zones are 
indicated by comers. They correspond to the critical Prandtl 
number Py^\\. In also varies with the nature of the waves. If 
A « 1, the Coriolis acceleration drives the dynamics and the 
perturbation generates quasi-inertial waves. If A » 1, it gener¬ 
ates quasi-gravity waves. This is the reason why we can observe 
two asymptotical curves on the plots. If we come back to the 
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Domain 


A ^ ^mn 

Pr » Pr,mn 

2 {nP + rP) E 

[nP + rP) E 

P f <$C P f^yyin 

rrP [rrP + rP)AEP-^ 

{irP +rP')EP-^ 

rP cos2 0 


Table 2. Scaling laws of the width at mid-height l„„ of the resonance 
associated to the doublet (m, n). is written as a function of the control 
parameters A, E and K = . 


case of inertial waves damped by viscous diffusion (jOgilvie & 
|Lin|[2004l l, 111 is proportional to E and does not varies with K. 
Thus, the peaks become thinner when E decreases that explains 
the behaviors observed in Figs. and In Appendix B, the 
reader will find color maps of the relative difference between 
the prediction of the analytical expression (Eq. [65] ), and the real 
theoretical result computed from the complete formula of 
(Eq. [3^ . They validate the expression of l\\ in asymptotical 
regimes and highlight transition zones where it cannot be used. 


Thanks to the scaling laws of Imn, it becomes possible to con¬ 
strain the variations of orbital parameters caused by a resonance 
of the internal dissipation. Eor instance, the amplitude Aa for the 


semi-major axis a of our two-body system (Auclair-Desrotour 


et al. 


2014 1 verify the law Aaja cx cx E in the regime 


ot inertial waves damped by viscous friction. 


3.3. Amplitude of resonances 

The height of a resonant peak also has an influence on the secular 
evolution of orbital parameters. Thus, we establish an analytical 
formula of the height of the peaks. Of course, the height of 
resonances depends on the tidal forcing f, which requires that we 
choose a form for it. Eollowing Ogilvie & Lin ( 2004| l, we take 
perturbation coefficients of the form: 


where we recognize the critical numbers and Pr.nm intro¬ 
duced in the previous section: 


^VISC _ 




i2 

nmn 


pI\a + -a„ 


m^rP (m2 -H rp) E{A+ Am„) {Pr -I- Pr,mn) 


(71) 


Thus, the height is characterized by the same asymptotic 
domains as the width. The frontiers between them depend on 
A„j„ and Pr,mn (Table [^. But it has different scaling laws. Eor 
example, when dissipation is led by viscosity, is inversely 
proportional to the Ekman number contrary to Imn- Therefore, 
the heights of the resonances increase if the Ekman number E 
decreases, which corroborates the plots of Eig. If dissipation 
is led by heat diffusion, then depend on both E and K. 
Of course, the scaling laws of are related to the spectral 
decomposition of the forcing. This latter partly determines the 
decay law of the heights with the wave number k: EPp^^'^ oc Ijk^. 
Here, oc Ijk^. But if f^,, did not depend on m and n, then 
the decay law would be oc l/k^. Consequently, with the 
forcing chosen for this study, the height of the resonant peaks 
becomes comparable with the non-resonant background rapidly. 
Therefore, few terms of the lowest orders are sufficient 
to plot the spectrum and there is no need to take into account 
a large number of harmonics to describe the dissipation properly. 


In the planet-satellite system studied by Auclair-Desrotour et al. 
(2014|l, we can thus write Aa/a with these new scaling laws: 


Aa 
— oc 
a 




(72) 


in the regime of inertial waves damped by viscous friction. 


_ , F 

fmn — l\ I Ti rjmn — 0? 1 ^ 1 ^^ — 0, (68) 

\m\ m 

and assuming {E, K] «c | ^/A, cos flj, we note that: 


P irUmn) ~ ^AnpK + (irP cos^ 6 + AirP^ isj , (69) 

so that we can compute simple expressions for the heights of the 
peaks of C™" and 


Domain 

A <3C Affifi 


P ^ P 

-* r -* r,mn 

AtiF'^E-^ 


irPrP {trP -t- 

m^rP {irP + tP)^ 

P P 

^ r -* nmn 

l6nFVcos'^0A-^E-^Pj 

^jiF^E-^P^ 

mP {nP + rP)^ 

nPrP {nP + tP)^ 


Table 3. Asymptotic behaviors of the height of the resonance of 
associated to the doublet (m, n). 


3.3.1. Dissipation by viscous friction 

Eirst, we look at the resonances of Naming 77",f their 
heights, we obtain the relation: 


tiy,. 


SttF^E cos^ 0 + ArrP'^ (rP cos^ 0 + AirP^ 

nPrP {nP + rP)^ [AnPK + {2rP cos^ 0 -i- AnP) Ef' 

(70) 


We plot the height of the main resonance (Eig.® as a func¬ 
tion of A, E and K, using the analytical scaling laws (Hq. [70|and 
Tab.®. The maximum of (Eig.[^ bottom-left panel) results 
frornthe inversion of the dependence on E, for Pr ~ Pr,n. So, it 
highlights the transition zone between viscous friction and ther¬ 
mal damping. As expected, the position of this maximum varies 
with A until A » 1, where Pr^w ^ 1 asymptotically. We recall 


in Appendix B). 


that Eq. ([70|) is relevant for small values of E and K (see Eig. B. 1 
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A = 10-4 
A = 10-4 
A= 10 , 
A =10-' 
A= 10, 
A= lo' 

A= 10^ 
A= 10^ 



A =10-4 
A = 10-4 
A = 10-f 
A =10-' 
A= 10, 
A= lO' 
A= 104 
A= 104 


log 10 E 


logio K 



A = 10-4 
A = 10-4 
A= 10 , 
A =10-' 
A= 10, 
A=10' 
A= lo; 
A = 104 



A =10-4 - 

A =10-4 . 

A= 10, 

A=10' . 

A = io; . 

A = 104 - 


log 10 E 


logio K 


Fig. 8. Width at mid-height In and height of the main resonance of as a function of the Ekman number E and its analogue for heat 
diffusion K for different values of A (in logarithmic scales). Top left: In - E. Top right: In - K. Bottom left: - E. Bottom right: - K. 


3.3.2. Dissipation by thermal diffusion 


An expression of 
deduced from Eq. ([3^, 


of the same form as Eq. (70i can be 


H' 


therm 


%nKF^ 


cos^ 6 + np-A 


rfi (m^ + tF-) \AnP-K + { 2 rP- cos^ 6 -H Airfi) E\ 


2 ’ 


(73) 


which can also be written as a function of the transition parame¬ 
ters Afyifi and Pr-^mni 


Domain 

A <^c! 

A ^ ^mn 

P ^ P 

1 p ^ 1 

ItiF’^E-^pY 

2,nF^A-^E-^p-Y 

rP {irP + rP-Y cos^ 6 

irPrP (irP + rpY 

P P 

1 p 1 p-^tfin 

^ttF^ cos^ 0A-^E-'^Pr 


nP {nP + rpY 

iii^rP (irP -1-11^)4 


Table 4. Asymptotic behaviors of the height of the resonance of 
^thcrm associated to the doublet (m, n). 


tl. 


therm 




Pr\A+ —A„ 


< 1 ^ (jfP + n^Y E(A+ A„„Y {Pr + Pr,mn) 


where H„„ ~ dissipation by heat diffusion, where 

(74) H„,„ ~ So, we introduce the ratio 


Thus, the asymptotic regimes of correspond to the four 

zones identified for l^n and (see Eig. |^. 


uvisc 

_ ^mn - p (A , A t 
“ utherm - 
^mn 


The transition zone is characterized by grnn 
3.3.3. Total dissipation where 


(75) 


1, OTl Pr 


pdiss 
^ r\mn'^ 


As and ^*enn j-jjg game resonance eigenfrequencies, 
their peaks are superposed. Therefore, when one of the dissi¬ 
pation mechanisms predominates at resonances, the observed 
peaks are all scaled by the law corresponding to this mech¬ 
anism. To establish the scaling laws of the total energy, we 
have to delimit the domains of dissipation by viscous friction. 


pdiss 
■* r;m« 


(0,A) 


1 

-A + Af^fi {&) 


(76) 


is the Prandtl number indicatiiw this transition zone. If Pr « 
Pf^„, then H„,„ ~ (Eig.g left panel); otherwise, H,„„ ~ 

Taking into account the new zones appearing with this 
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condition on Pr (Fig. right panel), we obtain the scaling laws 
for Hmn, the heights of resonances of the total dissipation given 
in Tab. For the sake of simplicity we introduce a new Prandtl 
number, 

(A) = max (A), (A)), (77) 


Assuming {E, K] « | VA, cos 0|, the previous expression be- 
comes 

(a -H COS^ 6 ») [ Cgrav (sll) A + Cin (sn) COS^ gl 

= 47rF2£e2 V-1 

® (A + cos^ef s] 2 +^(e,A,E,K) 

(82) 


that separates the regimes where Pr intervenes in scaling laws with 
from others. 


3.4. The non-resonant background 

The height of the non-resonant background i/bg is the magnitude 
of the energy dissipated inside the frequency range [winf, Wjup], 
in non-resonant areas. It only depends on the first term of the 
sum in the expression of ^ (Eq. [38] l. Indeed, each mode is 
related to a typical wavelength that increases when the wave 
numbers m and n decay. Thus, at the eigenfrequencies a> = comn 
{(m, n) e N^), the perturbation presents a characteristic pattern 
of wavelengths Ah - Elm in the x direction and A^, - Ljn 
in the z direction. But it is dominated by the lowest-order 
pattern m - n - I in the absence of resonance. If we had 
L ~ R, then our first mode would correspond to the large-scale 
hydrostatic adjustment of the flow in phase with the perturber, 
the equilibrium tide or non wave-like displacement introduced 
above (see e.g. |Remus et al.|2012a{|0gilvie|2013| l. 

In this framework, the height of the non-resonant background 
gives us information about the mean dissipation and the smooth 
component of the tidal quality factor (Q). Therefore, it plays 
an important role in the secular evolution of planetary sys¬ 
tems. Indeed, it is also necessary to compute the sharpness ratio 
E s //ii///bg intervening in the expression of Aa. Thus, in this 
subsection, we estimate by computing the term of the main 
resonance, at the frequency Wbg = (^ii + <r> 2 i) /2, that can 
be written: 


aihg - wii (1 + £ 12 ), 


(78) 


^ (0, A, E, K) = 20'^ (1 + 4ei2) - Aafijil + Isn) + 2y^. (83) 

The functions Cin (£ 12 ) and Cgrav (£ 12 ) are expressed by 


Cr(ei2) 

Cgt^(ei2) 


^2 -H 2ei2 -H £[ 2 ) (1 + £ 12 )^ 


(1 -1- £ 12 )“* 


(84) 


The coefficients p and y depend on E and K linearly (Eq.[5^. On 
the contrary, neither the Ekman number nor the dimensionless 
thermal diffusivity have an influence on cow and £ 12 . Therefore, 
the background linearly varies with the Ekman number if the 
condition 


If (0, A,E,K)\^[A-r cos^ of £^3 (85) 

is satisfied. To unravel what this implies on the physical param¬ 
eters, we need to study two asymptotic cases separately: 

1. A cos^ 0 (quasi-inertial waves), 

2. A » cos^ 0 (gravito-inertial waves). 


By deducing the main eigenfrequency a>u - \/04”+^o^^^)72 
from the dispersion relation (Eq. [25] l, we compute the relative 
distance between the two main resonances. 


with the relative distance between the two peaks. 


1 a)2i - wii 

£12 


2 a>ii 


V2- V5_ 

£12 ~ - 1 = - — Sir 


if A « cos^ 0, 


(79) 


£12 


2 V 2 - Vs 
V^ 


( 86 ) 


' £giav if A» cos 0. 


Note that if A = cos^ 0 (critical hyper-resonant case), the charac¬ 
teristic level of the non-resonant background is not defined and 
£12 = 0. Considering that the contributions of the main peaks are 
approximately the same, we write the characteristic height of the 
background: 


Hbg ~ 4f 11 (wbg). 


(80) 


3.4.1. Dissipation by viscous friction 

We first focus on the contribution of viscous friction. 


Thus, the relative distance between con and Wbg belongs to 
the interval [£in,£grav], £in and £grav being the boundaries cor¬ 
responding to the asymptotic cases A « cos^ 0 (i.e. the quasi- 
inertial waves) and A » cos^ 0 (i.e. the gravito-inertial waves) 
respectively. Numerically, £in ^ -0.183 and £grav ~ 0.132 and 
the asymptotic values of the coefficients and are: 


Qvisc ^ ^Visc ( ) ^ 32 87 

pvisc _ pvisc / \ _ ga 74 

'^grav ” '^grav \^*grav^ — 


(87) 


Erom Eq. we note that the dependence of the non-resonant 
background on E is linear only if: 


//b7^4f[r(cubg). 
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Domain 

A <SC ^mn 

A ^ ^mn 

1 y ^ L 

AnF^E-^ 


ni^rP- {m^ + rP-f 

nPrP {rrP + rpf 

P 

P r ^ P r\mn 

2nF^E-^Pf 

r r ^ r 

^nF^E-^P^ 

rP {nP + rpf cos^ 0 

irPrP {nP + rpf 

P r ^ Pr\mn 

^nF^ cos^ OA-^E-'^Pr 

rr ^ IP r,mn 

^jiF'^A-^E-^Pr 

niA {nP + rpf 

nPrP (nP + rpf 


Table 5. Asymptotic behaviors of the height of the resonance of ^ associated to the doublet (ra, n). Scaling laws correspond to the areas of 
Fig.|9] 


4 
3 
2 

^ 1 
Q. 

U) 0 

_o 

-2 
-3 
-4 

-4 -3 -2 -1 0 1 2 3 4 

log A 




-4-3-2-101234 

log A 


Fig. 9. Left: Zones of predominances for dissipative mechanisms. In the pink area (b), s> the dissipation is mainly due to thermal 
diffusion. In the white area (a), it is led by viscous friction. The transition zone corresponds to ~ P^'l\, where ~ f*'™. Right: Asymptotic 
domains with the predominance zones. Low Prandtl-number areas of Fig.|^are divided in sub-areas corresponding to the locally predominating 
dissipation mechanism. In a, b and f, ^ while in c, d and e, 


which corresponds to the quasi-adiabatic approximation met be¬ 
fore and assumed in the whole study. We thus finally obtain the 
expression of for inertial waves and gravity waves (Table 
|6ll. 


rrVlSC 

^bg 


= 47tF^E 


Cgl^A+q^'cos^g 

(A -H cos^ of' 


(89) 


Note that the background does not depend on the Prandtl 
number under the quasi-adiabatic approximation. Its level is 
only defined by the ratio A/cos^0 ~ A/An. So, the level of 
the non-resonant background is completely independent on the 
thermal diffusivity K in both asymptotic regimes. Besides, 

increases with E whereas decays. 


Note that thanks to our result for we now have a com- 

bg 

plete scaling law for the magnitude of the variation of the semi¬ 
major axis caused by a resonance of the dissipation spectra in the 
regime of inertial waves damped by viscous friction in a convec¬ 
tive planetary or stellar region ([Auclair-Desrotour et al.|[2()14[): 


Aa 

— cx /„ 






bg 


1/4 


oc E 


1/2 


(90) 


A All 

A » All 

A-Tia'^'^F^E 

in 

AnCf^F^A-^E 

cos^ 6 


Table 6. Asymptotic behaviors of the non-resonant background level 
of the spectrum. The case where A « cos^ 0 (A » cos^ 6) corre¬ 
sponds to quasi- (gravito-)inertial waves. 


It is interesting here to introduce a sharpness ratio H correspond¬ 
ing to the sensitivity of dissipation to the tidal frequency o). H is 
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A=10"l 
A=10 , 
A = 10 , 
A =10-' 
A= 10, 
A=10' 
A= lo; 
A= lO"* 



A = 10-^ - 

A = W-; . 

A =10"; . 

A _ 1A”1 . 


A= 10, 
A=10' 
A= lo; 
A= 10^ 


logic E 



A = 10-^ 
A = 10-^ 
A = 10"; 
A =10-' 
A= 10, 
A=10' 
A= 10^ 
A= lO"* 



W) 

o 



-9 -8 -7 -6 -5 -4 -3 -2 


A = 10-^ 
A=10 
A = 10-f 
A =10-' 
A= 10, 
A=10' 
A= lO, 
A= lO"* 


log 10 E 


logic K 


Fig. 10. Non-resonant background level and sharpness ratio as a function of E (with K = \Q ‘^) and K (with £ = 10 for different 
values of A (in logarithmic scales). Top left: - E. Top right: - K. Bottom left: - E. Bottom right: - K. 


defined as the ratio between the height of the resonance and the 
background level; 


analytical expressions of 
write the sharpness ratio; 


(Eq. 


70 1 and 

bg 


(Eq. 


82 1 , we 




(91) 


Therefore, the sharpness ratio estimates the relative magnitude 
of the variations of the dissipation due to the main resonant 
peak and compares the strengths of the dynamical and the 
equilibrium/non-wave like tides. A high sharpness ratio charac¬ 
terizes a peaked spectrum, and a low one a smooth dependence 
of the dissipation on u). As demonstrated in |Auclair-Desrotour] 
|et al.| p014| l, this parameter is important for the dynamical 
evolution of planet-moon and star-planet systems. Indeed, the 
variation of orbital parameters is correlated to S, the evolution 
of these laters being smooth in absence of resonances and erratic 
if E » 1 ( [Auclair-Desrotour et al.|201^ . 

We can define H for each of both contributions; 


^VISC 

T^visc _ 11 

uvisc 
bg 


^^herm 

and H*™ = 

^nerm 

'^bg 


(92) 


The transition zones of and are the same. Assuming 
the quasi-adiabatic approximation (Eq. [88|l, and considering the 


j ^2cos^-i-A^ ^Ah-cos^0^ Sj2 h-^( 0,A,£’,/f)j 
{AK + (2 COS^ H-A) [Cin COS^ 0 + Cgrav^] 


2^2 

*12 


(93) 


This expression may be simplihed in asymptotic domains (see 
Fig .0, where; 

x3 


"12 


(see Eq. 


\^i6,A,E,K)\ « (Ah-cos2 6») 

(ei 2 ) = (fiin) = (see Eq. |^, 

/^visc _ /^visc /„ \ _ ^vi 

'-"grav \^ll) — '^grav \^*grav^ — '-'gr 


85 I, 


"»V1SC 

grav • 


Thus, the sharpness ratio becomes; 


1 (l cos^ +a] (a + cos^ 

^Visc _ ^ _V_ /_\ _/_ 

'^[AK + (2 cos2 h-a) Ef [qq cos2 0 + qjfq]' 

Like the height of the resonances {ET^^), obviously 

presents symmetrical behaviors for gravito-inertial waves (Table 
0. It is inversely proportional to the square of the dominating di¬ 
mensionless diffusivity, E (for Pr » Pr,\\) oi' ^ (for Pr ^ Pr:\l)- 
So, the sensitivity of tidal dissipation to the frequency u in¬ 
creases quadratically when the diffusivities decrease (Eig. [T0| 
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bottom-left and bottom-right panels). Likewise, H increases with 
A quadratically when A » An. If A <Si An, then it is related to 
the colatitude 0. One will notice that Svisc does not depends on 
the forcing at all. Neither the amplitude F of the force, nor the 
form of the perturbation coefficients oc 1 / \m\ intervene 
in the expressions of Eq. (j9^ and ( |94| ). Given the link existing 
between the sharpness ratio and the tidal quality factor Q mod¬ 
eling dissipation discussed above, we can deduce the qualitative 
evolution of orbital parameters (Eq. [90l l. Indeed, if a planet is 
composed by a fluid with a large viscosity, the variation of the 
semi-major axis of its companion will be smooth. In the oppo¬ 
site case of a low viscosity planet, it becomes strongly erratic, 
with intermittent rapid changes of the orbit due to resonances 
(e.g Auclair-Desrotour et al.|2014ll. 


Domain 

A «c All 

A » An 

Pr » Pnll 

cos^ 6E^^ 

AE-^ 

Advise 

'^in 

o^visc 

■^'“'grav 

Pr « PrM 

cos® OA^^E^'^PI 

AE-^Pj 

/^visc 

^in 

o^visc 

grav 


Table 7. Asymptotic expressions of the sharpness ratio character¬ 
izing the spectrum of dissipation by viscous friction. 


3.4.2. Dissipation by thermal diffusion 


or by the interpolating function: 


Hi 


therm 

bg 


(A -H cos^ 0 )^ 


/ ^therm 
V'^grav 


A-hC*"™cos^ 


0 ). 


(100) 


Considering the asymptotic regimes, we obtain the scaling laws 
given in Tab. 


A An 

A » An 


47rC*"™£"A^"££;‘ 

EP? 

cos"^0 

Table 8. Asymptotic behaviors of the non-resonant background level 
77bg°*"' '■1'® dissipation by heat diffusion. The case where A <K An 

(A » An) corresponds to quasi-(gravito-)inertial waves. 


The sharpness ratio of dissipation by heat diffusion can be writ¬ 
ten: 


'-'therm 


A + -An \ en+^{e,A,E,K) 


4 £2 (A + An)" (Pr + Pniifd+enf’ 


( 101 ) 


Similarly as the height of the non-resonant background of expression which can be simplified in 
dissipation by heat diffusion is assumed to be approximatively: 


(95) 

Assuming the quasi-adiabatic approximation, we can write it: 


'-'therm 


Pi\A + -An 


^ EHA + An? {Pr + PrM? C*-” (sn) 


( 102 ) 


^ 4„F^EP-? 


{a + cos^ 0 ) (1 + £ 12 )" 

(A H- cos2 0)" s2 +^(g^A,E,K? 


C*“”(£i2) = (l + ^j , 


that will be a constant in asymptotical regimes, 

qherm ^ f^therm (g, j ^ 73 

C*r = C*-”(£g,av) = 73.10. 


Aj^p2pp-\ 

7/*"™ ~ -(ej2), 


bg 


(A + COs 2 0 ? 


in asymptotic regions where |log(A/Aii)| » 1, or by 


(96) 


where ^ has been defined above (see Eq. [83] ). We introduce the 
function C*"™ of ei 2 , 


therm 


7^^|A+-An 


4 £2 (A + An) (P, + £,;n)' + C^-^An) ’ 


(103) 


(97) 


if we use the asymptotic interpolation of the function C 


therm 


C*'=™(ei2(A)) 


(98) 


C therm a , ^therm a 
grav 21 + ^1 

A -I- An 


(104) 


When |log(A/An)| » 1, the height of the peaks can simply be 
expressed: 


(99) 


Using Eq. ( |103[ ), we establish the asymptotic scaling laws given 
in Tab. These laws are the same as the ones of Tab. which 
means that the variations of the energy with the tidal frequency 
do not depend on the dissipation mechanism in our local model. 
Therefore, the number of peaks, which is deduced from the 
sharpness ratio, will have the same scaling laws in both zones 
of Eig. |^(left panel). 
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Domain 

A «c All 

A » All 

Pc '3> P r ii 

cos2 6E^^ 

AE-^ 

g^therm 

■^^therm 

grav 

Pr « Pnll 

cos® OA^^E^^Pj 

AE-^Pl 

'^^therm 

'^in 

'^^therm 

grav 


Table 9. Asymptotic expressions of the sharpness ratio H*™ charac¬ 
terizing the spectrum of dissipation hy heat diffusion. 


3.4.3. Total dissipation 

The background level of denoted //bg, is the sum of and 
Qjygjj ^Yiat C*™ ~ C™" and C*™ ~ the ratio of 
the backgrounds is equivalent to the one of the heights. 


^VISC 

bg 

utherm 

bg 


^ 11 , 


(105) 


and we express scaling laws of i/bg in Tab. 10 for the regimes 
identified in Fig.[^(left panel). 

The asymptotic behaviors of the global sharpness ratio given in 
Tab. 11 are directly deduced from Tabs. andWe notice that 


^ - Mthe™ for all regimes. So, S ~ E' 

approximated by Tab. [7] 


and Tab. 


11 


can be 


Domain 

A «c All 

A » All 

Pr » 

' r;ll 

4;rC"^"£2£ 

47rq‘-T’2A-i£ 

COs2 0 

Pr P^^f, 

' r;!! 

47rC*™£2 

47rC*"™£’2A-2£f-;i 

cos^ 6 


Table 10. Asymptotic behaviors of the non-resonant background level 
i/bg of the total dissipation. The case where A <K An (A » An) corre¬ 
sponds to quasi-(gravito-)inertial waves. 


3.5. Number of resonant peaks 


In spite of the theoretical great number of resonances allowed by 
the condition given in Eq. ( [48] l (see also Fig. |^, the only effec¬ 
tive resonant peaks are the ones that exceed the level of the non¬ 
resonant background //bg (Figs. [3l and |4)l given by Eq. (|89]l and 
( |100| l (see also the discussion in |Kieutord & ValdettaroPoiu| l. 
This will tend to smooth the variations of the quality factor Q 
and the resulting evolution of the orbital parameters of a star- 
planet/planet-satellite system. Indeed, the resonances of highest 
orders are absorbed in the mean dissipation level of the resonant 
frequency range. The harmonics thus have to match a criterion 
to appear in the dissipation spectrum. This criterion simply cor¬ 
responds to the inequality; 


Hmn ^ ^bg» 

Article number, page 18 of 


24 


(106) 


where //„„ is given by Eqs. ( |7T] i and ( [7^ . The approximate wave 
number k introduced before allows to compute the typical order 
of the smallest peaks kc. As noted before, S ~ ~ 2*“=™, 

which implies that kc ~ ~ k*™. We replace the index m 

and n of the height //„„ by k and use Eq. (106 1 to obtain; 


^2 cos2 8 +A 


(a -H COs2 0^ 

^ el2+^iO,A,E,K)^ 

e],[AK + i2 

COs2 -hA) Ef' 

[Ci„ COs2 8 -t- CgravA] 


(107) 

This expression can be simplified in asymptotic cases. Assuming 
Eq. ([88]l, it becomes 


1 

[2 COs2 0 -1- a) 

(a -H COs 2 0j 

3 

1 

[AK +(2 cos^ h-A) Ef 

[C“cos2 0 + Cg“,A] 


(108) 

In Tab. we give the scaling laws of kc for the different 
asymptotic regimes. Taking into account resonances beyond 
this wave number does not change the global shape of the 
dissipation spectrum. In fact, in the situations corresponding 
to Figs |3 and|^ there is no need to go beyond k ~ 10 (Fig. 
|A.1| top-left and top-right panels). This is amply sufficient to 
model the dissipation realistically. The exponent 1/8 of kc in 
Eq. ( |107| i and ( 108| l is related to the spectral decomposition of 
the perturbing force. 


Domain 


A « A| 


P, » PrM 


P, « PrM 


1 

cos^ 6 \% 




cos° 6 \ 8 


C^A^E^P-^ 


A » All 




(_ 


Table 12. Asymptotic behaviors of the maximal order of noticeable 
resonances £• 


The wave number kc gives us the number of peaks (Eq. 61 1 
which is of great interest. As described above, we can assume 
Nk^ oc kc (Tab. [a and Fig.[^. Thus, we immediately deduce the 
scaling laws oiN^cc (Tab. [T^ from those of the rank of the har¬ 
monics (Tab. [T2|) in the four asymptotic regimes (see Fig.|^. In 
the same way as kc, Nk^ ~ ~ is given by the 

analytical expression; 


Nl 


^2 COs2 0 H- A 


(a + COs2 0^ 

%2^ -H^(0,A,£,/:)j 

[AK + (2 

cos2 -hA) Ef 

[Cin COs2 0 + CgravA] 


(109) 
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Domain 

A « All 

A » All 

Pr » 

' r;l 1 

cos^ OE^'^ 

AE-^ 

Advise 

'^in 

o^visc 

grav 

Pr « 

' r;l I 

Pr » F,,ii 

cos^ 6 E ^ 

Pr » P^'f, 

' r;ll 

AE-^P^, 

Q^therm 

'^in 

o^visc 

grav 

Pr Pr-n 

cos® 6 A-'^E-^Pl 

Pr P^'f, 

' r,\ 1 

AE-^Pj 

'^^therm 

'^in 

o^therm 

grav 


Table 11. Scaling laws of the sharpness ratio of f. They correspond to the areas of Fig.|^ 



A = 10-4 

A= 10 4 

A= 10 , 
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A= 10^ 
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A = 10-; 
A =10-' 
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A= 104 
A = 10- 


loglO E 


logio K 


Fig. 11. Number of resonances as a function of E (with A" = 10 **) and K (with £ = 10 **) for different values of A (in logarithmic scales). 

Left: - E. Right: - K. 


which is asymptotically equivalent to: 


1 

ft (2 cos^ 0 + a) (a + cos^ 0)^ 

^visc _ i-V-M- )_ - 

I 2 [AA: + (2 COs2 +A) Ef [C“ COs2 6 + C“avA] 

( 110 ) 

So, Nk^ oc £“02 fQj. inertial waves damped by viscous diffusion 
(see Fig. [m left and right panels). In this regime, the number of 
peaks decays when the Ekman number increases as can be ob¬ 
served in Fig.[^ The critical degree k^, the number of resonances 
Nk^. and the sharpness ratio E are related the ones to the others 
by the scaling equality: 


Nk,~kl~zK ( 111 ) 

in which, as pointed out above, the exponent 1/4 depends on 
the form of the coefficients of the perturbation, kc, Nk^, and H 
all measure the smoothness of the dissipation spectrum. As the 
change rate of orbital parameters is proportional to the energy 
dissipated inside the bodies, the evolution of an orbital system 
will be regular if they are small, erratic otherwise. 
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Pr « Frill Nt 


cos° 6 \4 


C°°A^E^P-^ 


Nk, 


■^'^grav^ r 


Table 13. Asymptotic behaviors of the number of peaks The same 

scaling laws are obtained for A*^. 


3.6. Super-adiabaticity (N^ < 0) 


The expression of (Eq. 38 1 allows us to compute dissipa¬ 
tion frequency-spectra if < 0 that corresponds to super- 
adiabaticity established by convection. Eig. shows how the 
spectrum of the energy dissipated by viscous friction evolves 
when the magnitude of an imaginary Brunt-Vaisala frequency 
increases. The graphs plotted correspond to pure inertial waves. 
We will not study here the particular scaling laws describing 
their properties, which are obviously not the same as the ones 
found for classical configurations. Nevertheless, we can observe 
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that strong values of |A| tend to flatten the spectrum. The height 
of the resonant peaks and the level of the background both decay. 
For exemple, the dissipation varies smoothly with the frequency 
when |A| = 10^ contrary to the case of positive N^. This situation 
will be studied in detail in a forthcoming work. 


4. Discussion 

Thanks to the analytical expression of the dissipated energy (Eq. 
[3^ , we have been able to compute the hydrodynamical scaling 
laws that link the properties of the dissipation to the dynamics 
of the fluid. They are summarized in Tab. [W and show that the 
related dependence of the dissipation onto the tidal frequency 
(a>) and the dimensionless characteristic numbers of the fluid 
(A, E and Pr) is tightly bound to the regime of tidal waves. In 
particular, they allow to explain the tendencies observed on the 
spectra presented by |Ogilvie & Lin| ( |2004| l. For example, it was 
noticed that the number of resonances grew if the Ekman number 
decreases in the regime of inertial waves dissipated by viscous 
friction (see Fig.[^. Now, we derive laws for our local model that 
demonstrate that it scales as Nk^ oc £“^2 Likewise, the height of 
the peaks has been analytically proved to be inversely propor¬ 
tional to E, contrary to the level of the non-resonant background 
for which we obtain //bg tx E. This explains why variations of 
the tidal frequency w strongly impact the dissipation when E 
tends to zero. Now, to understand how the regime of the waves 
determines the properties of the spectra, we have focused for in¬ 
stance on the sharpness ratio H. Table [l4| shows that this number 
always increases when the dominant dimensionless diffusivity, 
E or K, decreases. The difference between inertial waves and 
gravito-inertial waves resides in the linear dependence onto the 
stratification parameter A. This kind of comparison can be done 
for any other property similarly. In the light of the scaling laws, 
we identify the parameters that will have a real impact on the 
viscous dissipation and therefore, through the quality factor Q, 
on the long term evolution of planetary systems. 


5. Conclusion and perspectives 


Through a local approach, we have studied the physics of tidal 
gravito-inertial waves in stars and in fluid planetary regions. 
Their viscous and thermal dissipation is one of the main contrib¬ 
utor to the secular dynamics of spins and orbits. More precisely, 
it determines the so-called tidal quality factor Q (x) intervening 
in the dynamics ( Goldreich & Soter| 1966[ Auclair-Desrotour 
|et al.|[20T^ . This factor decreases when the dissipated tidal 
kinetic energy increases. Indeed, for an ideal elastic body, it 
takes very high values and the architecture of the orbital system 
is weakly affected by tides. At the opposite, if the body is the 
seat of a strong dissipation, the energy dissipated becomes 
important, the values of Q are smaller and systems evolve 
because of tides. This strong link between the secular evolution 
of a planetary system and tidal dissipation in their components 
therefore motivates studies to understand its dependence on 
their internal structure and dynamics. 


In the case of celestial bodies constituted of fluids, previous 
works highlighted a complex behavior (e.g. |Zahn|[T975i |Webb| 
|I980[[Ogilvie & Lin|2004[|2007) l. Indeed, liquid and gas layers 
do not dissipate energy like solid icy and rocky ones. The 
corresponding quality factor Q strongly varies as a function of 
forcing frequencies, fluid parameters (i.e. rotation, stratification 
and diffusivities), the geometry of the fluid container and bound¬ 
ary conditions and can vary over several orders of magnitude. 


Article number, page 20 of 


24 


This has important consequences for the evolution of spins and 
orbits (e.g. Witte & Savonije 1999[ Auclair-Desrotour et al. 


2014|l. In this study, the dissipation caused by viscous friction 


and thermal diffusion has been estimated in a local Cartesian 
section of a rotating body in order to understand and quantify 
these dependences in the whole domain of parameters relevant 
in stars and planets; in this framework, the Cartesian fluid box 
is complementary with global models (|0gilvie|[2005 Jouve & 


Ogilvie|2014[ Barker & Lithwick [2013 l2()14| l. The equations 


of momentum, mass conservation and heat transport are solved 
analytically with solutions expanded in periodic Fourier series. 
It yields the following conclusions: 


• Tidal dissipation resulting from viscous friction and thermal 
diffusion is highly resonant. Gravito-inertial waves are 
excited by tides at identified frequencies belonging to 
the interval [^inf^Tsup], where xm and Xsup are in a first 


approximation the inertial (2Q) and the Brunt-Vaisala (N) 
frequencies (for more precise boundaries, see Eq. and 
[Gerkema & Shriral|2005a] l. Thus, the fluid behaves like 
a bandpass filter. Inside the allowed frequency range, a 
batch of resonances can be excited by tides. Outside, the 
dissipation is weaker and varies smoothly with the forcing 
frequency that corresponds to the so-called equilibrium or 
non-wave like tide ( [Remus et al. 2012a| Ogilvie 2013) 1. In 
the resonant regime, a typical dissipation spectrum presents 
an organized structure of peaks which have different sizes 
and properties. These laters depend on the one hand on the 
spectral form of the forcing, in particular the amplitude 
of the harmonics which explains the fractal pattern of the 
studied case. On the other hand, they are narrowly bound to 
the fluid parameters: rotation (Q), stratification (N), thermal 
diffusivity (k) and viscosity (v). 


• There are four asymptotic regimes for tidal waves, charac¬ 
terized by the frequency ratio A and the Prandtl number P,-, 
and the positions of the transition borders are expressed as 
functions of the fluid parameters. The case when A « 1 
corresponds to quasi-inertial waves. They are mainly driven 
by the Coriolis acceleration even when > 0. If A » 1, 
then quasi-gravity waves are excited, the Archimedean force 
predominating on the Coriolis term. Besides, depending 
on the value of the Prandtl number Py, the dissipation is 
driven either by viscous friction (Py » Py-ii) or by thermal 
diffusion (Py « 

• For each of these asymptotic behaviors, all the properties of 
the dissipation spectrum can be expressed as functions of the 
fluid parameters: the positions, widths at mid-height, heights 
and number of resonances, the level of the non-resonant 
background and the sharpness ratio. These scaling laws, de¬ 
rived analytically from the model, can be used as zero-order 
approaches to constrain the resulting rotational and orbital 
dynamics. In this context, they constitute a tool to explore 
the physics of dissipation and offer a first basis to validate 
the results obtained from direct numerical simulations and 
global models. For instance, they can be used to compute 
scaling laws for the turbulent convective dissipation of tidal 
inertial waves in the convective envelope of stars and giant 
planets as a function of their rotation ([Mathis et al.|2014a)l. 


In a near future, the method applied in this work will be extended 
to other dissipative systems with for example differential rotation 
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Fig. 12. Spectra of in the case of pure inertial waves (£ = 10 and A" = 0) for A < 0. Top left: A = -0.1. Top right: A = -2. Bottom left: 
A = -10. Bottom right: A = -100. 


( Baruteau & Rieutord|2013| l or stratified convection ( [Leconte ^ 
|Chabrier||2012| l. Moreover, we will include magnetic field and 
take into account the corresponding Ohmic dissipation in addi¬ 
tion to viscous and thermal diffusions (e.g. [Barker & Lithwick] 
[2014[l. Tidal waves will b ecome magneto-gravito-inertial waves 
( Mathis & de Brye|201 l| l. This will introduce new fluid parame¬ 
ters (the Alfven frequency and the Elsasser and magnetic Prandtl 
numbers), new asymptotic behaviors and the corresponding scal¬ 
ing laws. Finally, non-linear interactions between tidal waves 
will be studied and characterized (e.g. [Galtier[[2003| [Sen et al.[ 
[2012] [Sen|20T3l [Barker & Lithwick|2013^ 
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Table 14. Scaling laws for the properties of the energy dissipated in the different asymptotic regimes. indicates the transition zone between 
a dissipation led by viscous friction and a dissipation led by heat diffusion. Ah indicafes fhe fransifion between tidal inertial and gravity waves. 
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Appendix A: Harmonic highest degree kc 

The expression of the highest harmonic degree, given in 
Eq. ( |108| l, is plotted as a function of the Ekman number and 
of the thermal diffusivity for various regimes from A 1 to 
A » 1. It corresponds to the upper bound of harmonic degree 
that is sufficient to represent all the resonances of the spectrum. 


Appendix B: Reiative differences rji and i]h 

The results provided by analytical expressions of Eqs. ( [65] l and 
fTT I have been compared to the ones given by the complete for¬ 
mula of the energy dissipated by viscous friction (Eq. [38] ). The 
color maps of Eig. |B.1| correspond to the computation of the 
width at mid-height In and height of the main resonance, 
summarized in Eig.[^ The relative differences rp and rjif are cal¬ 
culated using the expressions; 




Kana ^thl 


and rjH = 


uvisc _ uvisc 

l^ana ^th I 

uvisc 

^th 


(B.l) 


coming from Eqs. (651 and (iTT]), and 1* and 
being computed with Eq. ( [38| . The plots mghlight the asymp¬ 
totic domains, in dark blue, where the analytical formulae con¬ 
stitute a relevant approximation for the width and height of a 
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Fig. A.l. Rank of the highest harmonic degree as a function of E (with A = 10 '*) and K (with £" = 10 '*) for different values of A (in logarithmic 
scales). Left: “ E. Right: - K. 


resonance, and the critical transition zones where they cannot be 
applied. The light blue horizontal line corresponds to A ~ An, 
the hyper-resonant case. The color gradient in the regions near 
E ~ IQ-^ and K ~ lO^^ indicates the values of A, E and K for 
which the condition of the quasi-adiabatic approximation (Eq. 
[85] ) is not satisfied. The same can be done with the energy dissi¬ 
pated by thermal diffusion. 
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log E log K 

Fig. B.l. rj! and rju as a function of A, E and K in logarithm scales. White contours highlight critical zones where the quasi-adiabatic assumption 
({£, K] « { Va, cos nj) or the asymptotic condition (A t An) are not satisfied. Top left: Pi - (E, A). Top right: pi - (K, A). Bottom left: pn - (E, A). 

Bottom right: pn - {K,A). 
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